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ABSTRACT 

We use near ultraviolet and optical photometry to investigate the dust properties in 
the nearby starburst galaxy M82. By combining imaging from the Swift/UVOT in¬ 
strument and optical data from the Sloan Digital Sky Survey, we derive the extinction 
curve parameterized by the standard Ry factor, and the strength of the NUV 2175 A 
feature - quantified by a parameter B - out to projected galactocentric distances of 
4kpc. Our analysis is robust against possible degeneracies from the properties of the 
underlying stellar populations. Both B and Ry correlate with galactocentric distance, 
revealing a systematic trend of the dust properties. Our results confirm previous find¬ 
ings that dust in M82 is better fit by a Milky Way standard extinction curve (Hutton 
et al.), in contrast to a Calzetti law. We find a strong correlation between Ry and 
B , towards a stronger NUV bump in regions with higher Ry, possibly reflecting a 
distribution with larger dust grain sizes. The data we use were taken before SN2014J, 
and therefore can be used to probe the properties of the interstellar medium before 
the event. Our Ry values around the position of the supernova are significantly higher 
than recent measurements post-SN2014J (Ry mlA). This result is consistent with a 
significant change in the dust properties after the supernova event, either from disrup¬ 
tion of large grains or from the contribution by an intrinsic circumstellar component. 
Intrinsic variations among supernovae not accounted for could also give rise to this 
mismatch. 

Key words: galaxies: individual: M82 - galaxies: starburst - galaxies: stellar content 
- galaxies: ISM - ISM: dust, extinction - supernovae: individual: SN 2014J 


1 INTRODUCTION 

Extinction in astronomy describes how much light is ab¬ 
sorbed (or scattered) by dust grains in the medium be¬ 
tween the source and the observer. Locally, extinction is 
measured by comparing the light emitted from stars of sim¬ 
ilar spectral type, for which differential variations from in¬ 
terstellar dust can be assessed. Local measurements in the 
Milky Way reveal a general trend in the extinction law al¬ 


though with significant scatter (see, e.g. Fitzpatrick & Massa 
1990). The behaviour of the extinction law reveals the over¬ 


all properties of the dust grains in the interstellar medium 
I D raine||201l ). Closer to home, light scattered by the small 
air molecules in the atmosphere follow a typical oc A -4 
law that explains the blueness of the day-time sky. The 
typical Milky Way extinction law implies a shallower scal- 


dow, oc A 


Cardelli et al. 

1989 

Fitzpatrick 

1999 


effective law results from the integration of the contribu¬ 
tion from all grains in the dust component. A distribution 
with smaller grain sizes would make this law steeper, and 
a distribution with larger grains would induce a shallower 
dependence with wavelength, approaching the Mie regime 
(Draine 20111, where there is no wavelength dependence 


(i.e. clouds are white). One classical parameterisation of 
this smooth wavelength dependence is the ratio of total- 
to-selective extinction, Rv = Ay/E(B — V), where Ay is 
the extinction, measured in magnitudes, in the V band, and 
E(B — V) = Ab — Av is the colour excess, i.e. the additional 
reddening, measured as a B — V colour, caused by dust. In 
addition to this smooth trend in the extinction law, there 
are some features characteristic of a specific component. In 
our analysis - extended to the near UV and optical spectral 
window - the most re markable feature is the NUV bump at 
2,175A ( |Stecher|l969| , see Fig. I] Although the component 
responsible for this feature is not confirmed, it is possibly 
carbonaceous material, graphite, or PAH molecules (|Stecher| 
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& Donn 1965), although a more complex mixture is possible 
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Figure 1 . Top: Extinction law in the NUV (left) and NUV+optical (right) regions from [Conroy et al.| ( |2Q10] CSB10). Several cases of 
Ry and the NUV bump parameter B, are shown, as labelled. The CSB10 models on the left-hand panel all assume Ry = 3.1. In both 
figures, the dashed line shows for reference the Milky Way-standard extinction law (Fitzpatrick| 1999) F99). Bottom: In grey we show 
the response functions of the NUV filters of Swift/UVOT (UVW2, UVM2 and UVW1). The passbands of GALEX (FUV and NUV) are 
also shown on the left-hand panel, and the filter response curves of SDSS (u, g,r,i, z) are shown in the right-hand panel. For each filter 
the effective throughput is divided by the peak transmission. 


(|Bradley et al.|[2005|). NUV observations of the Milky Way 


extinction curve consistently find this feature (see, e.g. Fitz¬ 
patrick & Massa 1986 k In contrast, both lower metallicity 


regions, such as the Small Magellanic Cloud ( |Pei|19 92), and 
starbursting galaxies (Calzctti 20011, lack the NUV bump. 

An extinction law is a function that takes into account 
the wavelength-dependent obscuration due to absorption 
and scattering processes as light traverses a distribution of 
dust particles. One can therefore think of it as a foreground 
screen between the source and the observer. Such a scenario 
is a valid representation when observing stars in the Milky 
Way. However, when we extend the analysis of dust to dis¬ 
tant, mostly unresolved, galaxies, this process is complicated 
by the distribution of the dust, the variation of the compo¬ 
sition within the observed region, and the underlying stellar 
populations. It is common to refer to this integrated effect of 
dust in galaxies as attenuation. Galaxies with the same dust 
properties, hence with the same extinction, can yield differ¬ 
ent attenuation laws (see, e.g. Granato et al.[2000 1. For in¬ 
stance, the absence of the NUV bump in the attenuation law 
of starbursting galaxies (Calzetti 20011 can either be caused 
by an intrinsic lack of the component that causes this bump, 
or by the superposition of several stellar populations with 
different dust geometries ( |Panuzzo et al.|2007| ). Although it 
is beyond the scope of this paper to disentangle extinction 
and attenuation, one should bear in mind this subtlety in 
the interpretation of the observational constraints. 

This work focuses on the dust component in the nearby 
starburst galaxy M82 (NGC3034) located at a distance of 
3.3-3.5Mpc ( |Dalcanton et al. |2009|[Foley et al. |2014[) , with 


|son (20l2|). The ongoing strong star formation in this galaxy 
is believed to have been triggered by a close encounter with 
the more massive galaxy M81, located in the same group. 
UV light extends out from the central regions of the galaxy 
where star formation is the strongest, revealing the presence 
of dust entrained in the supernovas-driven wind. In |Hutton| 
feFal. (2014| hereafter HF14) we found that the dust obscu¬ 
ration in the NUV required the presence of an NUV bump, 
with the puzzling result that the prototypical starburst at¬ 
tenuation law of lCalzettl) ([2001) — which features no NUV 
bump — gave poorer fits than the standard Milky Way ex¬ 


tinction of Fitzpatrick (1999). This paper extends this issue 


by probing in more detail the dust properties of the inter¬ 
stellar medium (ISM) of M82, as a function of galactocentric 
distance. 

Literally, as the accepted HF14 paper went to press, an 
enthusiastic and diligently-led team of students at UCL’s 
University of London observatory discovered SN2014J, a 
type la supernova, within 1 kpc of the central starburst of 
M82 iFossey et al.|[2014). A number of papers followed the 


evolution of the supernova (see, e.g., Goobar et al. 2014 


Zheng et al. 2014 1 , and, most relevant to this paper, the 
properties of the dust (see, e.g. Amanullah et ak||2014 Fo 


a dynamical mass of ~ 10 1u Mq (Greco, Martini & Thomp 


|ley et aT~[ 1 2014 [ |Kawabata et ai.||2014[ |Patat et al. ||26i4| T 
We probe here the dust properties in the ISM of M82 pre- 
SN2014J, an issue of relevance to the understanding of the 
progenitor system of this supernova event. 

The structure of the paper is as follows: the NUV-to- 
Optical data used for the analysis are presented in (J2[ fol¬ 
lowed by the modelling of the dust attenuation law in 
The methodology followed to constrain the dust parameters 
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Table 1 . Parameter range of the grid of simple stellar population models used to derive the effective dust-related properties. 


Single Burst (1SSP): oc S(t 

- to) 


Observable 

Parameter 

Range 

Steps 

Age 

log(t 0 /Gyr) 

[-2,+0.9] 

32 

Metallicity 

log Z/Z© 

[-1.5,+0.3] 

8 

Total to selective extinction ratio 

Rv 

[0.5,4.5] 

32 

NUV Bump strength 

B 

[0,1.3] 

32 

Colour excess 

E(B-V) 

[0,1.5] 

32 


Number of models 8,388,608 


via stellar population synthesis is explained in © with the 
results shown in © along with comparisons with the recent 
constraints on the dust properties from SN2014J. Finally, © 
closes with the conclusions of our analysis. 


2 DATA 

We use the same dataset as in HF14. In a nutshell, deep 
NUV images were retrieved from the HEASARC archive of 
Swift/UVOT. The Ultraviolet/Optical Telescope (UVOT, 


Roming et al. 20051 is a 30 cm instrument working as a 


photon counter, with three passbands in the NUV region, 
straddling the position of the 2,175A bump. We refer the 


reader to Poole et al. (2008) and Breeveld et al. (2010) for 


details of the properties, mode of operation and data re¬ 
duction procedure of UVOT imaging. We retrieved imaging 
data of M82 with total exposure times of 10.2, 13.9 and 
8.7 ks in the UVW2, UVM2, and UVW1 passbands, respec¬ 
tively. The data were taken between 2008 and 2012, there¬ 
fore our photometry corresponds to the galaxy before the 
SN2014J event. In addition to the UVOT data, we retrieved 
optical images from the SDSS DR8 server (|Aihara et al. 


2011). All images were registered to a common grid with a 


0.5arcsec pixel size (i.e. the UVOT reference). The optical 
(SDSS) images were convolved to a common spatial resolu¬ 
tion with a Moffat profile. We note that the effective resolu¬ 
tion of the UVOT images is ~2.5 — 3arcsec. Photometry is 
performed within 5 arcsec (radius) apertures, following the 
standard calibration of the UVOT camera (jBreeveld et al. 


2010). The apertures were chosen to follow the disc of M82 


(see fig. 2 of HF14). We also retrieved the GALEX/FW 
images of M82 from the MAST archive, creating new sets of 
photometric measurements by convolving the UVOT data 
to the (poorer) resolution of GALE)© (~ 5 arcsec). There¬ 
fore, our combined dataset comprises photometry in the 
FUV, UVW2, UVM2, UVWl,u, g, r, i, z passbands. We re¬ 
fer the interested reader to HF14 for further details about 
the dataset. 

There are, however, two differences with respect to the 
HF14 data. Here, we make use of an updated estimate of the 
foreground (Milky Way) reddening towards M82. In HF14 
we took the value derived byjSchlafly & Finkb eliier| ([2011), 
as given in the NED database, E(B — V) = 0.140 mag. Fol¬ 


lowing the suggestion of Dalcanton et al. (20091, updated by 
|FoIey et al. ( |2014| ), we choose a lower estimate of 0.054 mag. 


1 Those lower resolution images were only used for the FUV — 
UVW2 colour 


This correction is due to the fact that M82 affects the read¬ 
ing of the large scale dust maps used for the determination 
of foreground reddening. This change affects the correction 
made to the observed photometry, and the derivation of the 
dust parameters, since in this paper, the amount of correc¬ 
tion from the Milky Way - which follows a Ry =3.1 [Fitz-| 
Patrick (19991 law - is significantly reduced. Nevertheless, 


we find that our results are compatible with HF14. In ad¬ 
dition, we use 3.3 Mpc for the distance to M82, for a more 


consistent comparison with recent papers (e.g. Foley et al. 


2014). This change only affects the translation from angular 


distances to projected physical distances. The photometric 
apertures adopted in this paper map into a projected 160 pc 
diameter. 


3 MODELLING THE ATTENUATION LAW 

We follow the parameterisation of the attenuation law of 
|Conroy et al.| (j20©j hereafter CSB10). This function de¬ 
pends on two basic parameters: Ry (defined in ©), and 
B, which controls the strength of the NUV bump. Fig. [l] 
shows a few realizations of the CSB10 function, along with 
the standard Milky Way attenuation law from |Fitzpatrick| 
(19991, assuming Ry = 3.1. Notice that within the NUV 


spectral window (left-hand panels), the UVOT filters sample 
the bump, with UVM2 located in the middle of the feature 
and UVW2/UVW1 straddling it on the blue and red sides, 
respectively. In addition, the GALEX/FUV filter includes 
an additional measurement blueward of the bump. The case 
B= 1 corresponds very nearly to the standard Milky Way ex¬ 
tinction curve, by construction. We follow a simple one-zone 
model for the dust component, assuming the light from the 
stellar populations is affected by a foreground dust screen. 


Other methods (e.g., Granato et al. 20001 would result in 


a more complex set of parameters, whereas we focus here 
on an “effective” attenuation law integrated to all the light 
within the aperture (see also Burgarcll a et al.|2005 |. 

In HF14, we compared the NUV and optical photome¬ 
try of M82 with a range of stellar population models, and 
two typical attenuation curves, the above mentioned Milky 


Way standard of Fitzpatrick (19991, and the law of Calzetti 


(20011 usually applied to star-bursting systems such as M82. 


We found that the latter - which lacks the presence of an 
NUV bump - gave poorer fits to the NUV photometry with 
respect to a Milky-Way law. In this paper, we explore in 
more detail this issue by use of the parameterization of 
CSB10. 

In addition to the NUV bump, the photometric data 
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Figure 2. Results from simulated data. The same method used for the analysis of the NUV+optical photometry of M82 is applied to 
a set of simulated spectra. Left: Direct comparison between input (i.e. true) and output (i.e. retrieved) parameters. Each row represents 
a set of 1,000 simulations where the adopted star formation history is a simple stellar population (SSP, top); an exponentially decaying 
star formation rate (EXP, middle) or a two-burst history (2SSP, bottom). The grey dots are individual results, whereas the big circles 
and error bars are median and RMS values, binned at fixed number of data points per bin. Right: The comparison is given with respect 
to the difference between retrieved and true values, expressed as a fraction of the uncertainty. The grey shaded regions and lines are 
the raw results, whereas the blue error bars correspond to a correction based on a calibration of the output using these simulations (see 
Tab. [ 2 ] and text for details). 



Retrieved Age (Gyr) Retrieved [M/H] 


Retrieved Age (Gyr) Retrieved [M/H] 


Retrieved Age (Gyr) Retrieved [M/H] 


Figure 3. Trend of the retrieved dust parameters with the age and metallicity derived from the same likelihood function. The offset is 
defined as An = 7r(in) — 7r(out), where n represents any of the three dust-related parameters. There is no significant bias in the dust 
parameters with respect to the underlying stellar populations. 


can be used to constrain the parameter Rv, which depends 
on the composition of the dust grains. A standard Milky 
Way law assumes Rv = 3.1, although with significant scat¬ 
ter, roughly between 2.2 and 5.8 ( Fitzpatrickf l999~I. Lower 
(higher) values of Ry will require a distribution of smaller 
(larger) dust grains, with respect to the Milky Way stan¬ 
dard. In HF14, we studied the extraplanar light originated 
from the reflection of the stellar light in the galaxy, that 


impinged on the dust grains above and below the galac¬ 
tic plane. These dust grains were entrained in the mate¬ 
rial ejected by the strong galactic winds triggered by the 
starburst. We fitted the photometric data assuming the ab¬ 
sorption/scattering processes are described by a power law 
(oc A -1 ), finding a value (x = 1.53 ± 0.17). Note, however, 
that this case pertains to the scattering of light from dust 
grains in the extraplanar wind region, as they are illumi- 
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OUTPUT - INPUT 


Figure 4. Histogram comparing the difference between the 
“true” (input) and the retrieved (output) dust-related parame¬ 
ters, for two choices of population synthesis models. The method¬ 
ology is the same in both cases, using a grid of SSP models from 
|Bruzual fe C hariot ( 2003| . The blue solid histograms show the pa¬ 
rameter retrieval accuracy when using a range of simulations con¬ 
structed from |Bruzual fe~Ch arlot (2003| models. The red dashed 
histograms correspond to simulations assembled from the mod¬ 
els of Marast.on] (2005[ ) ■ No significant bias is found (see text for 
details). 


nated by the galaxy. In contrast, we focus in this paper on 
the efective attenuation law in the main body of the galaxy. 

The right-hand panels of Fig. [l] show a wide 
NUV/optical spectral window, focusing on variations of Ry. 
Note how higher values of Rv - as expected in an ISM with 
larger dust grains - result in a flatter attenuation law, closer 
to the Mie regime. Therefore, the use of a wide spectral range 
(NUV and optical photometry) will allow us to constrain this 
parameter. 


4 CONSTRAINING THE DUST VALUES 


In order to constrain the dust attenuation parameters, we 
need to quantify the effect of the other, “nuisance”, pa¬ 
rameters related to the underlying stellar population. Our 
methodology involves a comparison of the observed NUV 
and optical photometry with a large grid of synthetic popu¬ 
lations from the models of Bruzual & Chariot (20031. Tab. [I] 
shows the grid of 8,388,608 models used in this paper. They 
correspond to simple stellar populations (SSP, representing 
a population with a single age and metallicity) for a stan¬ 


dard initial mass function (Chabrier 2003). Within our wave¬ 


length interval, dust simply introduces a smooth, long (spec¬ 
tral) range behaviour (aside from the NUV bump). There¬ 
fore it is possible to factor out the attenuation law from the 
different trends introduced in the photometry by the age 
and metallicity distribution of the populations. To prove this 


point, we ran three sets of simulations, each one with a differ¬ 
ent assumption about the star formation history: i) SSP (as 
above); ii) EXP: a composite of stellar populations follow¬ 
ing an exponentially decaying star formation rate, choosing 
the formation epoch, decay timescale and metallicity as free 
parameters of the populations; iii) 2SSP, a superposition of 
two SSPs, where the age of the young and old components, 
the fractional contribution, by mass, of the young compo¬ 
nent, and the metallicity are free parameters. These three 
choices encompass a wide range of both smooth and bursty 
star formation histories. This exercise serves to understand 
whether systematic biases are being introduced because of a 
methodology where the analysis is based on a grid of SSPs. 

The derivation of the model parameters follows a 
Bayesian approach, using a x 2_ based likelihood as a prob¬ 
ability distribution function (£ tx e~ Ax ^ 2 ) from which we 
derive the errors as the 68% confidence level of this distri¬ 
bution. Ax 2 is defined as x 2 — Xmmi an d X 2 is given by: 

T™OBS ™ MODELm2 

x 2 (m) - Y 

3 


J (7Ti)] 2 


a 2 (mj) 


where { 7 r;} represents all possible choices of the five pa¬ 
rameters listed in Tab. [l| rnf BS and j^N ODEIj are the 
observed and model magnitudes, respectively, as follows: 
{FUV, UVW2, UVM2, UVW1, u, g, r, i, z}. Finally, <r( mj ) 
is the observational errors for each magnitude. In addition 
to the parameters listed in Tab. [l] we include a “normal¬ 
ization” parameter that assumes the model SDSS g band 
magnitude lies somewhere between goBS ± 0.1 mag; where 
Sobs is the observed magnitude in the g band for each aper¬ 
ture. The extracted parameters are then computed using the 
likelihood as a probability distribution function, namely: 

f diV! ■ ■ ■ dn 5 £{n)Tri 

{7Vi) = fdm-dnsCin) (1) 

The priors applied in the analysis are simply the range of 
parameters explored, as shown in Tab. [T| (i.e. one can con¬ 
sider them as “flat priors”, giving equal probability within 
the range stated). 

Fig-i shows the results for 3x1,000 simulations cor¬ 
responding to each of the three formation histories. The 
simulated data explore a wide range of ages, metallicities, 
and dust properties, comparable to those found in M82 (see 
fig. 5 of HF14). Our analysis marginalizes over the stellar 
population properties, and produces the extinction-related 
constraints on - from left to right - Rv, B and E(B — V). 
The first two parameters are intrinsic to the effective dust 
attenuation law, whereas the latter can be considered as a 
normalization factor, reflecting the amount of dust present. 
The panels on the left compare input value (abscissa) ver¬ 
sus retrieved value (ordinate). An orange line in each panel 
gives the 1:1 correspondence. The grey dots are results from 
individual estimates, whereas the black dots and error bars 
mark the median and root-mean square values, respectively, 
when binned at fixed number of data points per bin. Note 
that no large bias is obtained, regardless of the use of differ¬ 
ent star formation histories (top to bottom). The derivation 
of the NUV bump strength only uses the UV photometry for 
the definition of the likelihood, since the optical photome¬ 
try is fully insensitive to the presence of this bump (Fig. [l] 
right-hand panel), and would therefore “wash-out” the con¬ 
straining power in the comparison. 
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Figure 5. Two-dimensional distribution of the retrieved values of the parameters that define the dust attenuation law used in this paper 
(see text for details). The ratio of total-to-selective attenuation (top-left), NUV bump strength ( top-right .) and colour excess ( bottom-left ) 
are shown as a colour-coded map, with the contiguous colour bars showing the mapping in each panel. The x-y coordinates are given as 
physical projected distances in kpc. The bottom-right panel shows a g-band image from SDSS-DR10 at the same scale. 


We use the comparison with the simulated data to de¬ 
rive a mean linear correction between input and output pa¬ 
rameters. The results of the fits are shown in Tab. [2] Since 
one cannot determine the best star formation history cor¬ 
responding to the real data, we decided to choose for the 
general correction the one that gave the best compromise 
among the three sets of star formation histories, i.e. we 
used the results from all three sets of simulations for the 
fit. The right-hand panels of Fig. [2] illustrate the validity of 
this correction and the retrieved error bars. They compare 
the difference between the “true” and the extracted value 
(“input” — “output” in the figure). This comparison allows 
us to determine whether the quoted error bars are valid. The 
shaded regions correspond to the original data, whereas the 
blue lines with error bars give the results for the corrected 
parameters. Note that the corrections are relatively minor 
with respect to the error bar, and no significant bias is ap¬ 
parent when the corrections are taken into account. Nev¬ 
ertheless, we will show in the next figures our results with 
and without this correction. The comparison with simulated 


data therefore confirms that it is possible to factor out the 
attenuation law from the properties of the underlying, unre¬ 
solved stellar populations, by use of a large library of simple 
stellar populations. 


We also need to assess a possible bias between the 
dust related parameters ( Rv,B,E(B — V)) and the under¬ 
lying stellar populations. In Fig. [3] we compare the bias 
An = 7r(in) — 7r(out) - where n represents any of the three 
dust-related parameters - as a function of the retrieved age 
and metallicity. Note that the input parameters are, by con¬ 
struction, uncorrelated. The figure shows that most of the 
comparisons give negligible dependence on age and metallic¬ 
ity, or on the functional form of the star formation history. 
This result is very important for the interpretation of the 
radial gradients of the effective attenuation law presented in 
the next section. We also note that the values we retrieve for 
the SSP ages are consistent with Hutton et al. (20141, where 
only two standard extinction laws were considered. These re¬ 
sults are also consistent with independent estimates, giving 
ages around 0.4-1 Gyr within the central 3 kpc, with older 
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Figure 6. Retrieved values for the ratio of total-to-selective ex¬ 
tinction Ry as a function of projected galactocentric distance. 
The grey dots represent individual measurements, with a charac¬ 
teristic error bar shown on the left-hand side of the figure. The 
red stars are measurements within 30 arcsec (i.e. 480 pc) of the 
position of SN2014J. The blue dots and error bars are the me¬ 
dian and root mean square values of the data, binned at a fixed 
number of data points per bin. The black line shows the same 
result if the retrieved parameters are not calibrated with respect 
to the simulations (see Tab. [2j. The orange cross and error bar 
correspond to the measurement in the region of SN2014J after 
the event (Amanullah ct al. 2014]) 


ages in the outer regions; and a colour excess E(B — V) ~ 0.4 
with a significant increase in the central regions (see, e.g. 


Mayya et al. 2006 Rodriguez-Merino et al. 2011). As re¬ 


gards to metallicity, we emphasize that broadband photo¬ 
metric analyses have large degeneracies with respect to this 
parameter. Our simulations in fact show that of all five pa¬ 
rameters given in Tab.[l] metallicity is the only one that fails 
at retrieving the input values, within a 0.5 dex uncertainty. 
However, we emphasize that such a result does not prevent 
us from deriving unbiased estimates of the other parameters 
(see Fig. [2j) , and consider age and metallicity simply as nui¬ 
sance parameters in the derivation of the dust-related prop¬ 
erties. In addition, we have done tests, imposing the typical 
metallicities expected in M82 (e.g. Origlia et al. 2004) as 
priors, and the extracted values of Ry, B, and E(B — V) 
are unaffected. A partial version of the parameters derived 
from this analysis are shown in Tab. [3} The complete version 
can be retrieved online. 

Finally, in order to assess the possible systematics 
caused by the specific population synthesis model used, we 
create an additional set of 1,000 simulations with the models 
of |Maraston] ( [2005] ). We note that these models apply widely 
different prescriptions, including a different treatment of the 
thermally-pulsing AGB phase, which can be prominent in 
populations of ages ~ 1 Gyr, therefore relevant for the study 
of M82. We emphasize here that the aim in this paper is not 


Figure 7. The observational constraints on the NUV bump 
strength, B , are plotted with respect to the projected galacto¬ 
centric distance. The symbols are the same as in Fig. [6] 


to constrain the age or metallicity composition of the pop¬ 
ulations, but to factor out the degeneracies caused by these 
parameters. In order to show that our results are robust 
against different prescriptions for the population synthesis 
models, we treat the simulated data - which use the |Maras-| 
ton (20051 - models, with the same grid of SSPs - which 
are based on the |Bruzual fe Chariot | ( |2003[ ) models. Fig. [4] 
shows the histograms of the difference between the “true” 
(i.e. the assumed value) parameter, and the retrieved one, in 
both cases (BC03: for the simulated data based on |Bruzual| 
fe Chariot] ( 2003]) models; CM05 for the data based on the 


Maraston (20051 models). We conclude that our derivation 


of the dust-related parameters is robust. 


5 THE ATTENUATION LAW OF M82 

Fig-0 shows the spatial distribution of the dust-related pa¬ 
rameters constrained with the photometry. The results are 
shown as colour-maps. The missing apertures in the distri¬ 
bution are masked out because of the presence of foreground 
stars. The bottom-right panel includes an SDSS-<? band im¬ 
age of the galaxy, shown at the same scale. The radial profile 
of the constraints on the dust properties of M82 are shown 
in Fig. [6] (Ry)', Fig. [7] (NUV bump strength); and Fig. [s] 
(colour excess). In all figures the grey dots correspond to 
individual measurements, with a typical error bar shown on 
the left-hand side. The median and root-mean-square val¬ 
ues of the data, binned at fixed number of data points per 
bin, are shown as black solid dots with error bars. Finally, 
the solid line without error bars show the result if the cal¬ 
ibration - presented in Tab. [2] - is not applied. Note the 
significant trend with galacto-centric distance. The central 
region of M82 - which due to projection effects include both 
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Projected Distance (kpc) 


Figure 8. The best fits to the colour excess, E(B — V), are shown 
against galactocentric distance. The symbols are the same as in 
Fig. [6] 


the starbursting region and the outer areas - has a more 
standard attenuation law, compatible wih the Ry = 3.1 and 
B = 1 values corresponding to the Fitzpatrick (1999) case. 
As we probe further out, the total-to-selective extinction ra¬ 
tio parameter, Rv, decreases, and so does the NUV bump 
strength, although we never find measurements compatible 
with a weak or non-existent bump. The colour excess also 
depends on galacto-centric distance, although in this case 
the correlation is more tightly linked to the variation in the 
metallicity and star formation distribution. 

We note that this result is robust against a possible 
bias caused by the inherent radial trend in age and metallic¬ 
ity found in galaxies. The previous section explored through 
simulations the possibility of a bias in the derivation of the 
parameters related to the dust composition, to find no sig¬ 
nificant systematic trends with age, metallicity or the dif¬ 
ferent functional form of the star formation history. There¬ 
fore, these trends are intrinsic radial variations of the effec¬ 
tive dust attenuation law ( Rv,B ), and the amount of dust 
(E(B - V )). 

Fig- E shows the correlations among the three dust- 
related parameters (using the same symbols as in the previ¬ 
ous figures). The range of values found within M82 is com¬ 
parable with the average estimates of Conroy et al. (2010) 
for a general sample of galaxies (Rv =2.0; B = 0.8). Our val¬ 
ues are also consistent with the sample of nearby galaxies of 


Burgarella et al. (2005), who find a range of values for the 


NUV bump strength, weaker than Milky Way standard, in 
the rangc|2]-B~0.5 — 0.6. Note the significant trend between 
Rv and B (top-left panel of Fig. ©), suggestive of a strong 


2 We translate their A(,u mp to the B parameterisation of Conroy 
|et al.| ( [2010 ) by a simple comparison of the effective attenuation 
laws. 


Table 2. Corrections to the data, derived from simulations (see 
fj4|. Col. 1 gives the parameter (tt) , cols. 2 and 3 are the slope and 
intercept of the correction; col.4 gives the scatter with respect to 
the full set of 3x1,000 simulated data points (given as an RMS); 
col. 5 is the bias, i.e. the offset between the fit and the data points. 


flTRUE = Q7r OBS + P 

7T 

a 

P 

RMS 

Rv 

1.030203 

-0.292036 

0.625 

B 

1.032269 

-0.111693 

0.225 

E(B - 

V) 1.027146 

-0.104647 

0.194 


link between the overall dust grain size distribution (larger 
grains causing higher values of Rv) and the carrier of the 
NUV bump strength. Although changes in the attenuation 
law can be produced by changes in the age-dependent ex¬ 
tinction of the geometry of the stellar/dust components, our 
study - based on a single galaxy - would suggest that most 
of the variations are caused by the intrinsic extinction law. 
A similar trend is found between B and colour reddening, al¬ 
though in this case the data mostly show that regions with 
high reddening ( E(B — V) > 0.5) have high NUV bump 
strengths (whereas less dusty regions can have a wide range 
of B, including B = 1). In a similar way, regions with high 
reddening are associated with higher values of Rv (top-right 
panel), corresponding to the presence of larger dust grains. 


5.1 The ISM around SN2014J in the 
pre-supernova phase 

Our data are especially relevant to explore the region around 
the recent SN2014J event. We use the position of the su¬ 
pernova from jTendulkar et a l. ( j2014j| to locate the rele¬ 
vant apertures in our dataset. In Fig. we show with 

red star symbols the estimates of Rv,B,E(B — V), respec¬ 
tively, measured in the apertures located within 30 arcsec of 
the SN2014J event, corresponding to a projected distance 
of 480 pc. Notice that those regions feature higher values 
of .Rv(~3-4) than the recent estimates determined post- 


SN2014J (Rv = 1.4 ± 0.1, Amanullah et al. 2014 orange 


cross). In general, photometric studies of the environment 
of type la supernova favour low values of Rv (Burns et 


al. 2014|), however the dispersion is quite large (Amanul- 
lah et al. 2015). For this particular case, even within the 


error bars and the dispersion of the data, it is not possi¬ 
ble to assume that post-SN2014J the dust properties in this 
region are compatible. Therefore, our results indicate that 
pre-SN2014J, this region had a dust content in the ISM sim¬ 
ilar to the Milky Way standard. Such a result implies the 
post-SN2014J observations should be explained either by a 
combination both from dust in the ISM surrounding the su¬ 
pernova and the intrinsic contribution from a circumstellar 


disk (Foley et al. 2014 although see Phillips et al. 20131, 


or by a significant change in the ISM towards the observed 
region. Alternatively, one could also expect that supernovae 
have intrinsic variations in colour that are not accounted 
for, or that the distribution of dust within the line of sight 
towards the supernova can produce the mismatch. 
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Table 3. Dust-related parameters and SSP ages derived from the analysis. Cols. 1 and 2 give the celestial coordinates of the central 
position of each aperture. Col. 3 is the corresponding projected galactocentric distance. Col. 4 is the total-to-selective attenuation ratio. 
Col. 5 is the NUV bump strength. Col. 6 is the colour excess, and col. 7 gives the base-10 logarithm of the SSP age, expressed in Gyr. 
Uncertainties are quoted at the 1 <t level. 


RA Dec. 

A R 

Rv 

B 

E(B - V) 

log Age 

J2000 

kpc 




Gyr 


148.95876 

69.67534 

0.365 

2.594 ± 0.089 

0.936 ±0.128 

1.025 ±0.042 

-1.979 ±0.060 

148.95539 

69.67785 

0.195 

2.612 ±0.095 

0.876 ±0.135 

1.061 ± 0.048 

-1.921 ± 0.153 

148.95203 

69.68039 

0.053 

3.626 ± 0.657 

0.793 ±0.147 

0.743 ±0.188 

-0.554 ±0.710 

148.94867 

69.68293 

0.171 

3.247 ± 0.570 

0.636 ±0.161 

0.513 ± 0.186 

-0.147 ± 0.603 

148.94530 

69.68547 

0.344 

2.777 ± 0.676 

0.537 ±0.150 

0.386 ± 0.208 

-0.133 ±0.703 


This is an excerpt from the central region. The full version of this table can be retrieved online. 


5.2 Are the dust properties of M82 different from 
other galaxies? 

We note that the trends presented in this paper are at odds 
with some recent results targeting our galaxy and M31. In 


Nataf et al. (20131 and Nishiyama et al. (20091, the dust at 


tenuation properties of the Milky Way are probed towards 
the Galactic bulge, finding a steeper law - i.e. lower Rv- 


In addition, Dong et al. (20141 find similar trends in the 


Andromeda Galaxy. We note that the populations towards 
the centre of M82 are significantly different from these two 
galaxies. The presence of a strong starburst will affect the 
properties of the ISM, and the dust. In addition, the nearly 
edge-on orientation of M82 along with its dusty ISM results 
in a complex interpretation of the results towards small pro¬ 
jected galactocentric radii. 

Another interesting issue is the strong correlation be¬ 
tween Ry and bump strength (top-left panel of Fig. [9]). A 
recent study of the NUV bump strength in high redshift 
galaxies (Kr iek fe Conroy| 26l3) suggests the opposite trend, 
namely a decreasing bump strength with increasing Rv- In 
contrast, detailed studies in our Galaxy do not seem to find 
a correlation (Fitzpatrick fe Massa||2007[ ), so there is still 
no clear indication whether such a correlation should be ex¬ 
pected in general. 

We emphasize here that our methodology is unlikely to 
produce a spurious correlation between these two param¬ 
eters. The simulations show that the (uncorrelated) dust 
parameters are retrieved without any significant bias or ar¬ 
tificial correlation among parameters. In addition, we note 


that our study mainly differs from Kriek & Conroy (20131 


in the sense that we explore internal variations within the 
same galaxy, rather than trends among galaxies. 


6 CONCLUSIONS 

We revisit the NUV and optical photometry of M82, pre¬ 
sented in HF14, to constrain in more detail the properties 
of the dust attenuation law. Through the analysis of sim¬ 
ulated data, we prove that it is possible to factor out the 
contribution from the underlying stellar populations, set¬ 
ting constraints on the NUV bump strength and the total- 



Figure 9. Correlations among the fitted dust-related parameters. 
The symbols are the same as in Fig. [fi] 


to-selective attenuation ratio Rv- We find a significant gra¬ 
dient with projected galactocentric distance (Fig. [6|, with 
the central regions showing a standard Milky Way attenu¬ 
ation law, progressing outwards towards an extinction with 
a weaker NUV bump and steeper Rv, more characteristic 
of a smaller grain size distribution, effectively tending to¬ 
wards the standard attenuation law of the SMC (Pei 1992). 
We should note here that our results correspond to a one- 
zone model, where the dust only contributes as a foreground 
screen. Changes in the geometry of the dust distribution can 
therefore cause variations on the dust parameters presented 
here (see, e.g., Witt & Gordon|[20001. 

As regards to the NUV bump strength, we find a sys¬ 
tematic decrease outwards. The weaker bump in the outer 
regions can be indicative of changes in the dust properties 
(Gordon, Calzetti, & Witt 1997|), although an age-dependent 
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attenuation law can produce similar effects (Panuzzo et al. 


2007). However, our observational constraints go in the op¬ 


posite way with respect to the models. A trend towards 
a weaker bump is found at large galacto-centric distances, 
away from the central starburst, where the dust content is 
lower, and where Rv is low, indicative of smaller dust grains 
(see Fig. [9|. Our results confirm HF14 in the sense that a 
Milky Way law, with a significant NUV bump is favoured 
with respect to a Calzetti (2001) law, in the central star- 
burst. 

Our work is most relevant to the analysis of the envi¬ 
ronment of the recent type la supernova SN2014J. Our data 
were taken before the event, and therefore probe the proper¬ 
ties of the dust of the ISM. At the position of SN2014J, we 
find a standard dust component, typical of the Milky Way, 
with 3 and B~l. Therefore, the recent measurements 
with lower values of Rv reflect either the presence of a cir- 
cumstellar component, or a significant change in the dust 
properties within the volume probed by the post-SN2014J 
measurements. 
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